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Reconstructing WIMP properties with neutrino detectors 



00 
O 

o 

\o 

(N 

S 1 

a; 

CD 

> 
OS 

o 

On 

m 

O 

t> 
o 



X 



Olga Mena 1 , Sergio Palomares-Ruiz 2 and Silvia Pascoli 2 
1 INFN Sez. di Roma, Dipartimento di Fisica, Universitd di Roma"La Sapienza" , P.le A. Moro, 5, 1-00185 Roma, Italy ana 
2 IPPP, Department of Physics, Durham University, Durham DH1 3LE, United Kingdom 

If the dark matter of the Universe is constituted by weakly interacting massive particles (WIMP), 
they would accumulate in the core of astrophysical objects as the Sun and annihilate into parti- 
cles of the Standard Model. High-energy neutrinos would be produced in the annihilations, both 
directly and via the subsequent decay of leptons, quarks and bosons. While Cerenkov neutrino 
detectors/telescopes can only count the number of neutrinos above some threshold energy, we study 
how, by exploiting their energy resolution, large magnetized iron calorimeter and, possibly, liquid 
argon and totally active scintillator detectors, planned for future long baseline neutrino experiments, 
have the capability of reconstructing the neutrino spectrum and might provide information on the 
dark matter properties. In particular, for a given value of the WIMP mass, we show that a future 
iron calorimeter could break the degeneracy between the WIMP-proton cross section and the an- 
nihilation branching ratios, present for Cerenkov detectors, and constrain their values with good 
accuracy. 

PACS numbers: 95.35.+d, 95.55.Vj, 14.80.-j, 95.85.Ry 



INTRODUCTION 

Cosmological and astrophysical observations provide 
growing evidence of the existence of Dark Matter (DM) in 
our Universe. One of the favored candidates is a weakly 
interacting massive particle (WIMP), a new stable neu- 
tral particle with typical mass in the 10 GeV-TeV range 
which interacts weakly [l| . These particles arise naturally 
in various extensions of the Standard Model of Particle 
Interactions (SM), e.g. the lightest neutralino in super- 
symmetric models with R-parity conservation 1], stable 
scalars in little Higgs models [2j, the lightest Kaluza- 
Klein excitation of the SM fields in models with univer- 
sal extra-dimensions [3[ . Direct searches for WIMPs look 
for the recoil energy of nuclei due to interactions with 
the DM particles passing through the Earth Q. WIMP 
signatures can also be found indirectly, namely by ob- 
serving gamma-rays, positrons, anti-protons and neutri- 
nos produced in dark matter annihilations [Ij. As the 
annihilation rate scales with the square of the dark mat- 
ter density, indirect signals are expected to be stronger 
where DM accumulates, as, for instance, the galactic cen- 
ter and astrophysical bodies. High-energy neutrinos are 
produced in the annihilation either directly or by sub- 
sequent decay of SM fermions and bosons [f| and their 
spectrum could provide information on the WIMP mass, 
the WIMP-proton cross section and the branching ratios 
of the various annihilation channels, and hence, on the 
DM properties. 

Searches for high energy neutrinos can be performed in 
Cerenkov neutrino detectors/ telescopes laul, as Super- 
Kamiokande (SK), AMANDA, IceCube, but no positive 



evidence has been found so far. However, these detec- 
tors are unable to fully reconstruct the neutrino energy 
and therefore cannot provide information on the neutrino 
spectrum. In the case of small detectors with large photo- 
multiplier coverage, like SK, the size limits the maxi- 
mum energy for which the events are contained, so only 
through-going muons are of interest for energies of tens of 
GeV or higher. In the case of neutrino telescopes like Ice- 
Cube, least granular, the size limits the energy threshold 
for detection. In both cases and for the energies of inter- 
est, these detectors can only count neutrinos above their 
energy threshold. On the other hand, different concep- 
tual designs for detectors for future neutrino long base- 
line experiments show that Magnetized Iron Calorimeter 
DetectorsJMIND) H, or its variants as MONOLITH [| 
or INO [lfjL L iquid Argon Time Projection Chambers 
(LArTP C) fill and Totally Active Scintillator Detectors 
(TASD) [l2| are able to measure the energy and direction 
of the initial neutrino with good precision at energies of 
tens of GeVs. 

In the present letter, we restrict our analysis to MIND, 
which are conceptually similar to the existing MINOS de- 
tector 13( but with a mass one order of magnitude larger. 



Due to its ability for a precise muon charge discrimina- 
tion, they have been proposed as ideal detectors for fu- 
ture neutrino factories. For the same reason, these detec- 
tors have also been considered as atmospheric neutrino 
observatories, opening up the possibility of determining 
the neutrino mass hierarchy and the small mixing angle 
#13, due to matter effects in the oscillations of upward- 
going atmospheric neutrinos [14j. Here, we reveal a novel 
aspect of MIND |l5[ , showing that they might constitute 
excellent observatories to detect neutrinos from WIMP 
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annihilations. For a given value of the WIMP mass, 
we study the neutrino signal from WIMP annihilations 
in these detectors, and exploit their energy and angular 
resolution in order to reconstruct the neutrino spectrum 
and, ultimately, to study DM properties [16j |. By mea- 
suring the branching ratios of annihilation into different 
channels and the WIMP-proton cross section, these mea- 
surements could provide a unique opportunity for getting 
information on the WIMP couplings with the particles of 
the SM. 



NEUTRINOS FROM WIMPS ANNIHILATIONS 

When a WIMP goes through a celestial body, as the 
Sun, it might interact elastically with the nuclei and 
get scattered to a velocity smaller than the escape ve- 
locity, remaining gravitationally trapped in the body. 
It will then undergo additional scatterings, settling in 
the Sun core, giving raise to an isothermal distribution. 
The WIMPs so accumulated can annihilate with other 
WIMPs into SM particles, as quarks, leptons and, if kinc- 
matically allowed, gauge and Higgs bosons. Typically in 
the Sun, for sufficiently high capture rate and annihila- 
tion cross section, equilibrium is reached and the anni- 
hilation rate r ann is related to the capture rate Cq as 

r„ 



and depends on the WIMP mass and on the branching 
ratios into the various channels: 



I/2C0, which, for the Sun, is given by [Tj 
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where p\ oca \ is the local WIMPs density, wi oca i is the ve- 
locity dispersion of WIMPs in the halo, a is the WIMP- 
proton cross section, and m DM is the WIMP mass. The 
capture rate can be different by several orders of mag- 
nitude for spin-dependent and spin-independent interac- 
tions. Typically, for most neutralino models as well as for 
Kaluza-Klein DM, the spin-dependent cross-section dom- 
inates over the spin-independent one. The opposite hap- 
pens for scalar WIMPs. These scattering cross sections 
are constrained in direct DM searches [18j, which put 
very strong bounds on the spin-independent ones. Since, 
in the Earth, the abundance of nuclei with odd mass 
number is extremely small, the neutrino flux from WIMP 
annihilations is strongly constrained by the bounds on 
the spin-independent cross section and would be too low 
to be interesting for detection in MIND, LArTPC or 
TASD. Thus, for our purposes we will focus on WIMP an- 
nihilations in the Sun, for which the much more weakly 
constrained spin-dependent cross section may play the 
dominant role. 

High-energy neutrinos would be produced directly or 
indirectly via the decay of other products of the annihila- 
tions [l| . A broad spectrum of neutrinos is so generated 
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where the sum includes the possible annihilation channels 
with spectrum dNi/dE v and branching ratio BR^, and 
R is the Sun-Earth distance. If the DM is a Majorana 
fcrmion, the annihilation amplitude into fermion pairs is 
proportional to the fermion mass, so that the neutrino 
channel is irrelevant and annihilations are dominated by 
heavy fermions, 6s, rs, cs and, if kinematically allowed, 
is. Neutralinos can annihilate into fermions as well as 
gauge bosons and Higgs and the dominant channels are 
controlled by the neutralino composition Typically, 
for neutralinos lighter than the W boson, the annihila- 
tion into bb pairs gives the main contribution, with harder 
neutrinos from rf. For higher masses, the gauge bosons 
channels are allowed and can have the highest branching 
ratio. Kaluza-Klcin WIMPs can annihilate directly into 
neutrinos |3[ with a branching ratio of few per cent and 
would provide a specific signature with a peak in the neu- 
trino spectrum at E v = rn DM . Their main annihilation 
channels are charged leptons and light quarks with the 
neutrino flux generated dominantly by the rf mode. In 
the following, we will treat the branching ratios as free 
parameters as their exact values depend on the type of 
particle considered and on its couplings. 

Annihilations into t's and gauge bosons produce a hard 
spectrum which peaks around 0.4 m DM (0.5 m DM ) for v &lil 
(y T ) for the former channel and even at higher energies 
for the latter [l9j]. Quarks hadronize producing a large 
number of mesons and baryons which subsequently de- 
cay producing a softer neutrino spectrum. For instance, 
neutrinos from 6s and cs have a spectrum which peaks at 
energies ~ 0.15 m DM or lower, depending on the WIMP 
mass, and drops to zero around E v ~ (0.6-0.7) m DM [ljj ]. 
Light quarks hadronize mostly in pions which, as muons, 
get stopped before decaying and therefore produce only 
low energy neutrinos. Once produced, neutrinos propa- 
gate in the Sun, being absorbed via charged-current in- 
teractions and loosing energy due to neutral-current ones. 
These effects, in addition to neutrino oscillations between 
different flavors, need to be taken into account (lij . 



RECONSTRUCTING DM PROPERTIES 

High-energy neutrinos can be detected in neutrino de- 
tectors and stringent bounds have already been set 
and will be further improved by neutrino telescopes 
However, as these experiments cannot reconstruct the 
neutrino energy, they can only provide limited informa- 
tion on the neutrino spectrum. In order to reconstruct 
the neutrino spectrum, the analysis presented here ex- 
ploits the energy and angular resolution of MIND, which 
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allow not only to measure precisely the energy and direc- 
tion of muons produced by the charged current interac- 
tions (CC) of v ,. an d but possibly also the energy of 
the electrons [8lf20l|. produced by v e [v e ) interactions. In 
addition, these detectors will be able to reconstruct the 
energy and angle of the hadron shower, so that a good 
energy and direction resolution for the incoming neutrino 
will be achieved. We show that the degeneracy among 
different parameters can be broken by using the muon 
and electron energy spectra. In particular, a larger an- 
nihilation rate r ann could always be traded by a smaller 
branching ratio to a hard channel (t + t~) in absence of 
information on the spectrum. 

In the following we will consider a light WIMP can- 
didate with m DM < 80 GeV and therefore, we will focus 
on its dominant annihilation modes for these masses, i.e., 
t + t~ (hard) and bb (soft) channels. Note that the con- 
tribution from cc is similar to the one of bb but with a 
lower normalization [l9j |. This technique might not be 
competitive with neutrino observatories (AMANDA and 
Icecube) for WIMP masses ^ 100 GeV, since the typ- 
ical sizes of future MIND, LArTPC and TASD might 
not allow a precise measurement of high lepton energies. 
However, it should be noted that, even for high WIMP 
masses, these detectors could be useful for detecting the 
low energy tail of the neutrino spectrum (23| . In our 
study we assume the mass of the WIMP to be known and, 
for definiteness, we take m DM = 50 GeV and 70 GeV. 

Let us now comment on our assumption of taking a 
given value for the WIMP mass. By the time the detec- 
tors discussed in the present article are available, infor- 
mation on the WIMP mass might be already available 
from the LHC (Large Hadron Collider) and, possibly, a 
future ILC (International Linear Collider) as well as di- 
rect and indirect dark matter searches. With the start of 
LHC in 2008, the neutral WIMP candidate for the dark 
matter can be indirectly detected in an event looking for 
missing energy and missing transverse momentum. By 
a detailed determination of the kinematics of the quarks 
and leptons, typically it might be possible to measure the 
WIMP mass with a 10% accuracy (see, e.g., Ref. [24]). 
The prospected error depends strongly on the model of 
physics beyond the SM and the values of the parameters. 
For example, for the SPSla benchmark point, the LSP 
mass could be determined to be m DM = (96 ± 5) GeV at 
the LHC and m DM = (96 ± 0.05) GeV at the ILC 0]. 
For less favorable points in the parameter space, worse 
accuracies will be reachable. The value for the mass de- 
termined in collider searches could be used as input in our 
analysis, assuming that the observed WIMP corresponds 
to the dark matter particle. 

A model-independent determination of the WIMP 
dark matter mass could also be obtained in direct DM 
searches which look for the nuclear recoil in an inter- 
action between a WIMP and a nucleus in the detec- 



if m DM < 100 GeV . For instance, a superCDMS-like 
experiment, with exposures of 3 x 10 3 (3 x 10 4 ) [3 x 
10 5 ] kg days and a spin-independent cross section close 
to the present bound ersi = 10~ 7 pb, would be able 
to measure a light WIMP mass with an accuracy of 
- 25% (15%) [2.5%] [26|]. Similar conclusions can be 
drawn for the case of spin-dependent cross sections with 
comparable statistics. Nevertheless, these measurements 
are affected by the uncertainties in the WIMP veloc- 
ity distribution and in the local WIMP distribution. A 
method which is independent of the WIMP density and 
of the WIMP-nucleon cross section has been also pro- 
posed [271 ]. By comparing the recoil energy spectrum in 
direct searches which use different detector materials, it 
is found that a WIMP mass smaller than 50 GeV can be 
determined with a 1-tr error of 20% if a total statistics of 



tor 25| . The recoil spectrum is strongly mass-dependent 



500 events is available [271 ]. 

Information on the WIMP mass could also be found 
from indirect dark matter searches. Some information 
on the WIMP mass could also be obtained from the 
angular distribution of neutrino-induced muons in large 
water- Cerenkov detectors (28[. On the other hand, as the 
WIMP density in the center of galaxies is predicted to be 
enhanced, a high rate of WIMP annihilation are expected 
and produce sizable fluxes of SM particles, as neutrinos, 
photons, quarks and leptons. High energy photons can be 
produced either directly or via the hadronic interactions 
in the interstellar medium. In the first case, the gamma 
spectrum presents a very clear signature with a peak at 
the WIMP mass, allowing its precise determination. In 
the case of gamma-rays from primary hadrons, the spec- 
trum has a sharp energy cut-off in correspondence with 
m DM . It has been estimated that GLAST could constrain 
the WIMP mass at the 25% level, for m DM = 100 GeV, 
and future Imaging Atmospheric Cerenkov Telescopes, 
with sufficiently small energy threshold, could detect the 
peak from direct annihilation into photons [29| . 

In addition, the energy information in the WIMP neu- 
trino signal in the detectors considered in the present ar- 
ticle could allow, in principle, to infer directly the WIMP 
mass by exploiting the cutoff in the spectra. Since the 
typical sizes of future MIND, LArTPC and TASD might 
not allow to have sufficient statistics at the tail of the 
spectrum, this measurement could be very challenging. 
A detailed analysis needs to be performed, but it is be- 
yond the scope of the present study Q . 

We have computed the number of electron and muon 
neutrino- induced CC events for m DM = 50 and 70 GeV, 
by using the evolved neutrino spectra from Ref. f]~9j ] . We 
have, conservatively, considered nine energy bins (5 GeV 
bins from an energy of 5 GeV for m DM — 50 GeV and 
7 GeV bins from an energy of 7 GeV for m DM = 70 GeV) , 
where the low energy cut is adopted to avoid the at- 
mospheric neutrino background. With good statistics, 
the precision in the measurement of m DM , for m DM > 
50 GeV, depends on the muon energy resolution of MIND 
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at energies E v ~ 50 — 100 GeV, expected to be ~ 5 — 8% 
at 50 GeV and slightly worse for higher energies [2(| ■ 

For an ideal detector with perfect lepton angular res- 
olution, the atmospheric neutrino background would be 
reduced to the fraction of events within the angular size 
of the Sun (~ 6.7 x 10~ 5 sr). However, the lepton does 
not travel in the same direction as the incident neutrino 
except when averaged over many events. We conserva- 
tively take the root mean square (rms) spread in direction 
between the incident neutrino and the lepton (in radians) 
as 



1 GeV 



E„ 



(3) 



where E v is the energy of the incoming neutrino. There- 
fore, a very conservative approach is to sum over the 
atmospheric neutrino background contributions within a 
region of opening angle rms . 

The expected number of muon (electron) neutrino- 
induced events in the i-th bin is computed as a func- 
tion of the WIMP-proton scattering cross section and of 
the branching ratios into b + b~ (soft channel) and t + t~ 
(hard channel), 
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where A is the energy bin width, 9 is the angle which 
measures the relative position of the Sun with respect 
to the detector location, Nt is the number of available 
targets, Vdot is the total volume of the detector, t is the 
exposure time, <fi is the evolved (anti)neutrino spectra, 
<7 CC is the CC (anti)neutrino cross section and V^, e is 
the volume of detector available for the neutrino to in- 
teract. For electron-like events, which are all contained, 
this fiducial volume V e is simply the detector volume Vdet ■ 
For muon-like contained-events, it depends on the detec- 
tor geometry and on the muon range in iron R U (E U ). For 
a detector with cylindrical shape, it is given by [2l| 



V^E^ff) = 2hr 2 arcsin [ \l 1 - R ^ E ^ sm 2 ■ 



4H 
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where r = 13 m and h — 20 m are the detector radius and 
height we have considered. We approximate the muon en- 
ergy by its averaged value, which for the energy range of 
interest is (E^) = 0.48 -Ej, [22[. Also notice that the vol- 



ume of detector available is a function of the azimuthal 
angle 9, which varies from to 7r in 12 hours. Thus, for 
the assumed geometry, the effective volume will be maxi- 
mal when the Sun is at the horizon. On the other hand, if 



E 
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a good energy resolution is achieved up to some energies 
even for partially-contained events (by means of using 
the bending of the muon trajectory due to the magnetic 
field), the fiducial volume to be considered is Vdot- We 
will consider both cases, when partially-contained muon- 
like events are not included in the analysis and when all 
muon-like events with the interaction vertex in the de- 
tector contribute. 

The expected number of muon and electron neutrino- 
induced events is then fitted performing a x 2 analysis. 
We define 

(6) 

where C is the covariance matrix and only statistical er- 
rors, being the dominant ones, have been considered. 
Ni tCC is given by Eq. and n^ a represents the simu- 
lated muon (electron) data taking into account the to- 
tal muon or electron atmospheric neutrino background 
within the i-th energy bin and integrated over a solid 
angle Afi = 27r# rms , and is given by 

n i>a = Smear(7V iiQ + N i>atmoa ) - A^ atmos , (7) 

where the function Smear indicates that Gaussian or 
Poisson smearing has been applied following the Monte 
Carlo techniques of Ref. [30] to mimic the statistical un- 
certainty. 

In Figs. Q] and [21 we illustrate the simultaneous ex- 
traction of the scattering cross section of the WIMP off 
a proton and the branching ratio into t + t~ assuming a 
50 kton MIND with the geometry described above. In 
Fig. Q] we consider the case of a WIMP mass of m DM = 
50 GeV and the input simulated values are er = 7xl0~ 3 
pb and BR r + r - = 20%. We have assumed 15 years of 
data taking and have only included electron-like and con- 
tained muon-like events. The dotted-dashed red lines 
depict the 90% confidence level (CL) contours assuming 
that no energy information is available: as previously 
anticipated, a larger annihilation rate could always be 
mimicked by a lower branching ratio into t + t~ . Notice 
that there will be a continuous region of degenerate solu- 
tions, implying the impossibility of extracting either the 
different branching ratios or the WIMP-proton scattering 
cross section. In this scenario, the measurement of the 
WIMP mass is of course out of reach. Instead, if we now 
include the information on the neutrino energy spectrum, 
we obtain the 90% CL contours depicted by the solid blue 
line. The strong correlation between a and BR T + T - (or 
BRb+f,-) is broken. The reason is simple: the differential 
neutrino spectrum for each possible annihilation channel 
(hard or soft) has a characteristic shape, different from 
one channel to another. For this case, m DM = 50 GeV, the 
spectrum for bb peaks at ~ 7 GeV and represents roughly 
20% of that of t + t~, peaked at ~ 20 GeV. We also il- 
lustrate the results for the ideal case for which the atmo- 
spheric neutrino background only contributes within the 
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FIG. 1: The solid blue (dotted-dashed red) contours denote the 
90% CL limits, for 2 degrees of freedom, for the simultaneous ex- 
traction of the WIMP-proton cross section and the annihilation 
branching ratio into the t + t~ channel with (without) energy in- 
formation by using only fully-contained events with an exposure of 
750 kton ■ yrs (see text for details) . The simulated nature values are 
m DM = 50 GeV, a = 7x 10~ 3 pb and BR T+T _ = 20%. The dashed 
magenta lines (inner contours) illustrate the case (for 90% CL) in 
which the atmospheric neutrino background is integrated only over 
the real size of the Sun. We also show the 90% CL excluded region 
by SK 0. 
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FIG. 2: The solid blue (dotted-dashed red) contours denote the 
90% CL limits, for 2 degrees of freedom, for the simultaneous ex- 
traction of the WIMP-proton cross section and the annihilation 
branching ratio into the t + t — channel with (without) energy in- 
formation by using only fully-contained events with an exposure 
of 750 kton ■ yrs (see text for details). The simulated nature val- 
ues are m DM = 70 GeV, a = 5 X 10" 3 pb and BR T+T _ = 10%. 
The dashed magenta lines (inner contours) illustrate the case (for 
90% CL) in which partially-contained events are also included in 
the analysis. We also show the 90% CL excluded region by SK Q. 



real angular size of the Sun (dashed magenta lines). As 
can be seen, by considering a more detailed and realis- 
tic expression for the angular resolution than that given 
by Eq. ([3]), our results would not be significantly im- 
proved. In Fig. [3 we show the 90% CL contours for 15 
years for m DM = 70 GeV with the input simulated values 
a = 5 x pb and BR r+r - = 10%. The solid blue 

line represents the case when partially-contained muon- 
like events are not considered, whereas the inner dotted 
magenta line depicts the case when we also include this 
type of events to perform the analysis. In this figure, we 
have considered the atmospheric neutrino background in- 
tegrated over a half-cone aperture given by Eq. ^ . It is 
important to note that by adding the partially-contained 
events, in order to achieve the same results as for the 
case of only contained events, the total number of years 
could be reduced from 15 to 10 years of data taking, the 
minimum running time of these detectors in the context 
of neutrino factories. In both figures, the equivalent 90% 
CL SK limit @ is included. 

As both figures show, for these relatively low WIMP 
masses and large spin-dependent cross sections, the de- 
termination of the annihilation branching ratios and 



WIMP-proton cross section could be achieved. Notice 
that a positive signal of neutrinos from WIMP annihi- 
lations in the Sun in these detectors will immediately 
point to a spin-dependent cross section. In LArTPC and 
TASD similar results could be obtained. In fact, we have 
checked that considering only the or the v e signal does 
not modify the results substantially, a part from a small 
change due to the smaller atmospheric neutrino back- 
ground for v e . A detailed analysis for different WIMP 
masses and type of detectors is in preparation [23| . 



CONCLUSIONS 

Searches of high-energy neutrinos from WIMP anni- 
hilations in the Sun could constitute a powerful probe 
of WIMP properties. While the total neutrino flux is 
controlled by the annihilation rate (proportional to the 
WIMP-proton cross section), the shape of the neutrino 
spectrum depends on the WIMP mass and on the branch- 
ing ratios of different channels. Ccrcnkov neutrino detec- 
tors/telescopes are counting experiments and can only 
provide limited information. On the contrary, we have 



6 



shown that, by exploiting the good energy and angu- 
lar resolution at energies of tens of GeV, future MIND, 
LArTPC and TASD could have the capability to recon- 
struct the neutrino spectrum, and could provide impor- 
tant information on the WIMP mass, its annihilation 
branching ratios and WIMP-proton cross section. In 
particular, our analysis shows that, for a given value of 
the WIMP mass, the degeneracy between the WIMP- 
proton cross section and the branching ratios into soft 
and hard channels could be broken and important infor- 
mation on the WIMP dark matter couplings might be 
obtained. This information should be combined with the 
constraints and measurements from direct and indirect 
DM searches to probe the nature and properties of DM 
particles and to test that the WIMP candidate found in 
collider searches does indeed constitute the observed dark 
matter of the Universe. 
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